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Abstract 



Effects of thermal conductivity on the dynamics of first-order phase transitions 
are studied. Important consequences of a diff'erence of the isothermal and adiabatic 
^ ' spinodal regions are discussed. We demonstrate that in hydro dynamical calculations 

fl^ i at non-zero thermal conductivity, k ^ 0, onset of the spinodal instability occurs, 

when the system trajectory crosses the isothermal spinodal line. Only for k = it 
occurs at a cross of the adiabatic spinodal line. Therefore ideal hydrodynamics is 
not suited for an appropriate description of first order phase transitions. 
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Q . 1 Introduction 

o 



There are numerous phenomena, where first-order phase transitions occur be- 
K> ! tween phases of different densities. In nuclear physics various first-order phase 



H [ transitions may take place in the Early Universe [I], in heavy ion collisions [2|3 

and in neutron stars (e.g., pion, kaon condensation and deconfinment or chiral 
phase transition [HE])- At low collision energies the nuclear gas-liquid (nGL) 
first-order phase transition is possible [2|6|7f8|9] . The experimental status of 
the latter is discussed in Refs. p^ll|12j . At high collision energies the hadron 
- quark gluon plasma (hQGP) first-order phase transition may occur (see e.g. 
Ref. [3]). Within the hydrodynamical approach dynamical aspects of the nGL 
and hQGP transitions have been recently studied in Refs. [T5|14|15] . An im- 
portant role of effects of non-ideal hydrodynamics has been emphasized in 
Ref. [T3fTi] . However in these works density fluctuations were carefully ana- 
lyzed, whereas thermal transport effects were only roughly estimated. 

Dynamical trajectories of the system in heavy-ion collisions can be approxi- 
mately determined by constant initial values of the entropy per baryon (i.e. 
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entropy density per net baryon density). This conclusion follows from the anal- 
ysis of heavy ion collisions at energies less than several GeV/A performed in 
an expanding fireball model (e.g., see Ref. [ISj), from ideal hydrodynamical 
calculations in a broad energy range (see e.g. Refs. PTPS] ) and also from 
transport calculations (e.g., see Ref. [T9]). 

It is well known, see Refs. [T5|20|21 ]. that, at least in the mean field approx- 
imation, the isothermal spinodal (ITS) line and adiabatic spinodal (AS) line 
are different. From the exact thermodynamic relation one obtains 
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Here T is the temperature, P is the pressure, n is the density of the con- 
serving charge (here the baryon density), m is the baryon mass, s = s/n is 
the entropy per baryon, cy is the specific heat density at the fixed volume V , 
u1 = m~^{dP/dn)s and m^ = mr^{dP/dn)T- The value Ug has the meaning of 
the adiabatic sound velocity (at s =const). Namely, this quantity character- 
izes propagation of sound waves in ideal hydrodynamics. As we show below, 
in non-ideal hydrodynamics at a finite value of the thermal conductivity, k, 
the propagation of sound waves is defined by the interplay between ut (at 
T =const) and Ug- Conditions m^ = and m^ = define T(n)-curves: the 
ITS line (at fixed temperature T) and AS line (at fixed s), respectively. The 
maximum temperature points on these lines are the critical temperature T^r 
in case of T=const and the adiabatic maximum temperature Tp^max in case of 
s =const. In the mean field approximation, the specific heat density cy has 
finite non-negative values. Thus one has m|. < on the AS line, where Ug = 0. 

Calculations performed in mean field models show that Tcr > Tp^max- In 
Refs. [7|9] and numerous subsequent works, crossing of the AS line by effective 
trajectories s=const was considered as the starting point for the multifragmen- 
tation following an adiabatic expansion of the system. In contrast, in Ref. [8] 
and later Refs. pT|22f23] it was assumed that the multifragmentation may 
appear already at higher temperatures, when the ITS curve is crossed by the 
s=const lines. The argumentation of the former group of works can be traced 
back to the fact that in ideal hydrodynamics instabilities indeed arise at the 
AS boundary (see extensive discussions in Ref. [15]). Isothermal instabilities 
at first order phase transitions were studied in details in Refs. [T^fT^ . 

In mean field theories, difference between values of Tcr and Tp^^nax is usually 
rather large. For example, within the NJL model one obtains Tp^^nax ~ ^cr/2 — 
45 MeV (see Ref. [21]) for the chiral restoration transition. In the model used 
in Ref. [15] Tp^max = 3X^^/5 ~ 100 MeV for hQGP phase transition. Note 
that simultaneously with the low value Tp^^nax ~ (50 ^ 100) MeV a sufficiently 
high value of npmax should be reached. This looks unlikely to be achieved in 



heavy-ion collisions. Hence, if instability occurs for T < Tp^max, one could face 
difficulties in exploring the hQGP spinodal instabilities in actual experiments 
with heavy ions. 



It was argued in Ref. [21] that above mentioned difference of Tcr and Tp^rnax is 
probably an artifact of the mean field approximation and both values of the 
critical temperature should coincide provided thermal fiuctuation effects are 
included. Indeed, it is known (e.g., see Ref. (2^), that inclusion of fluctuations 
may result in a weak divergence of the speciflc heat. For example, for the three- 
dimensional Ising universality class cy ~ [{T — Tc)/Tc]~", where a ^ 0.12. It 
can be concluded from Eq. ([1]), that for divergent cy both temperatures T^^ 
and Tp^^^ should coincide. Although this observation is fully true, it seems 
legitimate to rise some questions. 

• First, thermal fluctuation effects are strong enough only in a critical region, 
where the mean fleld theory breaks down and system properties are governed 
by non-trivial critical exponents. For many systems the critical region is a 
rather narrow region near the critical point, while contrary examples are 
known as well [21] . Although it is not known a priori how large is the hQGP 
critical region, there is an indication in favor of the hQGP critical region 
being small [25]. Only if the system trajectory passes trough the critical 
region, one can hope that fluctuation effects appearing during evolution of 
the system may lead to T^^ ~ Tp^^^. 

• Second, even if the system trajectory passes through the critical region and 
T^^ coincide with Tp^^^, the ITS and the AS lines may not collapse into 
one line even within the critical region because the fluctuation contribution 
to the speciflc heat diverges only at the critical point and then sharply 
decreases. Thus we should understand crossing of which of two lines (the 
ITS or AS line) results in the spinodal instability. 

• Third, the divergence of cy at the critical point is due to long-scale fluctu- 
ations, which need a long time to develop. Therefore the time, which the 
system spends in a vicinity of the critical point in the course of the heavy 
ion collision, might be insufficient for the developing of such fluctuations, 
see discussions in Refs. [T3lll4ll26ll271. 



Moreover, if the initial collision energy is such that the system trajectory 
passes outside the critical region, i.e. entirely in the mean fleld regime, the 
thermal fluctuation effects do not contribute at all. Then the question arises 
how one can treat the problem purely within the mean fleld approximation. 
This is actually an important question, since most of equations of state used 
in the description of heavy ion collisions disregard the fluctuation effects. 

In Refs. [T3pl] we focused on the study of dynamical effects of the viscosity 
and the surface tension. The aim of the given paper is to consider the influence 
of thermal conductivity effects on the dynamics of flrst-order phase transitions 
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Fig. 1. Phase diagram of the Van der Waals equation of state, T(n)-plane. Bold 
sohd, dashed and dash-dotted curves demonstrate boundaries of the Maxweh con- 
struction, the spinodal region at T =const and s =const, respectively. The short 
dash lines show adiabatic trajectories of the system evolution, see explanation in 
the text. 




Fig. 2. Phase diagram for the Van der Waals equation of state in P{n) plane, see 
explanation in the text. 

and to answer above risen questions. Below we find out that in non-ideal hy- 
drodynamics (with non-zero thermal conductivity, k ^ 0) instabilities show up 
at the crossing of the ITS line and only for k = 0, e.g. in ideal hydrodynamics, 
they appear at the AS hne. 



2 General setup 



To simplify further considerations we describe the dynamics of a first-order 
phase transition by means of the standard system of equations of non-ideal 



non-relativistic hydrodynainica ^ I: the Navier-Stokes equation, the continuity 
equation, and the general equation for the heat transport: 



mn [dtUi + (uV)Mi] = -ViP (2) 

2 



+Va 



ri 



V...V,.-i...vuUa.d.vu 



dtu + div(nu) = 0, (3) 



T 



ds 

- + div(.u) 



: div(fi:VT) (4) 

/ 2 \2 

+r] iVkUi + ViUk - -5jfcdivuj + C(divu)^. 

Here r] and C, are the shear and bulk viscosities; u is the velocity of the fluid 
element; s is the entropy density; k is the thermal conductivity; d is the di- 
mensionality of space. The short-scale fluctuations can be taken into account 
by adding noise terms on the right sides of Eqs. ([2111]). In what follows we 
neglect these contributions. 

Before we proceed further we aware the reader, that Eqs. ([2l|4]) are derived in 
the first gradient order approximation, see Ref. |28]. In the kinetic regime the 
typical inverse relaxation time of the system t~2 and the typical inverse spatial 
scale fckin characterizing the collective processes fulfil the inequahties r~/ <^ ech 
and fckin ^ Pch, see ^29j, where ech and pch characterize microscopic processes. 
For quasi-equilibrium nucleon system with T > ep (ep is the nucleon Fermi 
energy) one has ech ~ T and Pch ^ Pt = v2rnT. The hydrodynamic regime is 
reached if the typical time scale and the typical length scale characterizing the 
collective motion are larger than the corresponding kinetic values, Xhydro ^ T"rei 
and /hydro ^ 1/^kin- In what follows we assume that these inequalities are 
fulfilled. 

An approximate validity of the adiabatic expansion (i.e. conservation of spe- 
cific entropy density, s, along the fluid flow) in heavy ion collisions [19] means 
that (i) K, is rather small, or the system is spatially homogeneous and (ii) i] 
and ( are small, or the velocity of the expansion of the system is very low 
(with the guaranty for absence of shock waves). Note that in the general case 
solutions of Ref. [13] , which we will exploit below, are obtained at a non-zero 



^ A short survey of some results generalized for relativistic case is done in Appendix 
B. 



viscosity therefore we will not additionally assume smallness of the viscosity 
coefficients, rather we exploit smallness of the velocity u. 

The best known example to illustrate principal features of a first order phase 
transition in the mean field approximation is the Van der Waals fiuid. Expres- 
sions for thermodynamical quantities in this particular case are presented in 
Appendix A. Trajectories for the expansion of the uniform system, obtained at 
the assumption s ~ const, are shown in Fig. 1 on the plot T /T^r = fin/ricr)- 
The upper convex curve (MC, bold solid line) is the boundary of the Maxwell 
construction, the next bold dashed line (ITS) shows the boundary of the ITS 
region and the lower bold dash-dotted curve (AS), the boundary of the AS 
region. The short dashed lines crossing the regions under these three curves 
are calculated at two values of the entropy per baryon: Scr and Sm- Different 
scenarios of the dynamics for first-order phase transitions in expanding sys- 
tem at s ~ const have been considered previously in Ref . [S] . The super-cooled 
vapor (SV) and the overheated liquid (OL) regions are between the MC (bold 
solid curve) and the ITS (bold dash) curves, on the left and on the right re- 
spectively. For Scr > S > 5mc25 where Scr corresponds to the critical point and 
the line with smc2 passes through the point n/ucr = 3 at T = 0, the system 
traverses the OL state (region OL in Fig. 1) and then the ITS region (below 
ITS line) and the AS region (below AS line). For s > Scr the system trajectory 
passes through the SV state (the region SV in Fig. 1) and then the ITS region. 

In Fig. 2 we show the same regions as in Fig.l, but in the plane P/Pcr ~ n/ucr 
for P > 0. The line Sm corresponds to the maximum pressure at the AS 
curve. The critical point is determined by {dP/dn)T = 0, {d?'P/dn'^)T = 0. 
With the help of Eq. f l37|) of Appendix A from these conditions we find Ucr = 
1/(36), Tcr = 8a/{27b). The maximal pressure point at AS line is defined by 
{dP/dn)s = 0, {d'^P/din?)s = 0. Thus using Eq. (HUj) we obtain: np^rnax = 
1/(96) = Ucr/^i Tp^jnax = 128a/(12156) = 16Tcr/45, see Fig. 2. The maximal 
temperature point on the AS curve corresponds to Tmax = 3/5Tcr, Umax = ^cr, 
see Fig. 1. Another relevant point on each line s = const is the turning point, 
where E{n, s = const) has a minimum as a function of the density, n. In this 
point the pressure P{n,s = const) equals zero and the fireball expansion is 
slowing down. The lines labeled Sn and Sj2 cross respectively the AS and ITS 
lines at P = 0. Note, that with realistic equation of state, P = might be 
reached at the nGL first-order phase transition, see Ref. [20j, but not at the 
hQGP transition. 

We stress that we use the Van der Waals equation of state only for the illus- 
trative purpose. Our main conclusions do not depend on the used equation of 
state provided there is a first order phase transition. 

We may expand all quantities near an arbitrary reference point (nr,Ti.), or 
(rZr, Sr). Analytical solutions can be found, provided the system is in the vicin- 



ity of the critical point {ricr, Tcr), or near the maximum point {np^max, Sp^max)- 
Therefore considering the problem in (n, T) variables it is convenient to take 
(ucryTcr) as the reference point and considering the problem in (n, s) vari- 
ables, to take ijip^rnaxiSp^rnax) as the reference point. The following remark is 
in order. As we have mentioned, in this paper we exploit the mean field ap- 
proximation. Thus using {ncr,Tcr) as the reference point we actually assume 
that the fluctuation region [ricr =F ^n^, Tcr =F ST^) is very narrow and we are 
able to chose the reference point outside this region (but nevertheless very 
close to the critical point). Also, mentioned problem is avoided provided one 
considers the system at the time scale, being shorter than that necessary for 
a development of long-scale critical fluctuations. 

We construct the generating functional in {6n = n — ricr, ST = T — Tcr) or 
in {6n = n — np^rnax, Ss = s — sp^rnax) Variables, the Landau free energy, such 
that 5iSFL)/5i5n) = P - Pf + P^c- 



r d^x j cm[V{Sn)]'^ \m'^{5n)'^ Xv'^m{6n)'^ 1 

^~y^j 2 ^ 4 2 ^^J' ^^ 

where e = Pf — Pmc is expressed through the pressure on the Maxwell con- 
struction. The maximum of the quantity e is e™ = AXv^ / {3\/3) , therefore we 
may count e from this value supposing further that e = ±7^6'" with < 7^ < 1. 
The first term in Eq. (^ is due to the surface tension. It produces a surface 
pressure Psurf (see Eq. (TTUj) below). 



In case of the Van der Waals form of the equation of state, see Appendix A, 
expanding the pressure of the uniform matter in Sn = n—ncr, ST = T~Tcr, for 
\6n\/ncr ^ 1, \6T\/Tcr <^ 1, i.e. in the vicinity of the critical point {ricr, Tcr), 
we obtain 



SP[n,T] = fo 



9ST5n 9Te,((5n)3 3 ,^ 



(6) 



The last explicitly presented term in Eq. ([6]) is unimportant if T(?^) ~ const, 
since the addition of any constant does not change equation of motion ([2]). 

Comparing Eq. ([6]) with the value, which follows from Eq. ([5]), we find relations 
between coefficients: 



v\T) = -46Tnl.myZ 



cr. 



'-r = tI^. (7) 



Expanding the pressure with respect to 6n = n — np^rnax, Ss = s — sp^rnax in 



the vicinity of the maximuin of P{n) at fixed s we obtain 



6P[n, S] = —rTp^max 



U{5nf 114 

OSdn + — ^ h —^np^jnaxOS + 



Comparing this expression with the value, which follows from Eq. (|5]), we find 
relations 



v\s) = 

'^P,max 






P,max 



m^ 



(9) 



Further, following Refs. [I3II1], we use smallness of the velocity u(t, r) of the 
growth/damping of the density and temperature fluctuations. We linearize 
all terms in hydrodjTiamical equations in the velocity u , the entropy density 
6s = s — Sr and the temperature 5T = T — T^ variables. Here index "r" stays 
for an arbitrary reference point. 



Excluding a dependence on the velocity, u, from Eqs. ([2]), (|3]) we obtain p^l4|30] 



md^Sn = A 



P + Psurf -Pf + n;^ (dT], + Cr) dt6n 



(10) 



where Psurf is the surface pressure; d = 2{d — l)/d, d = 4/3 for spherical 
geometry, 6n = n — n^. We count P from its value in the final equilibrium 
state (reached at t — t- oo). Namely this difference, P — Pf, has the meaning of 
the thermodynamical force driving the system to the final equilibrium state. 

Eq. ( TTU]) should be supplemented by the equation for the heat transport (|1]), 
which after linearization owing to Eq. ([3|) simplifies as 



dtSs — Sr{nj.) ^dtSn = Kj.A5T. 



;ii) 



The variation of the temperature is related to the variation of the entropy 
density s[n, T] by 



6T ~ T,{cv,r)~^ {Ss - {ds/dn)T,M) , 



(12) 



where, as in Eq. ([T]), cy^j. is the specific heat density. 



3 Instabilities in spinodal region 



3.1 Effect of finite thermal conductivity 



Since in this section tlie "r" -reference point can be taken arbitrary, we sup- 
press index "r". To find solutions of hydrodynamical equations we present, cf. 
[TMillTS] . 



6n = 6nQexp [yt + ipr] , 

6 s = 6soexp['yt + ipr], 

T = T> + (5Toexp[7t + ipr]. 



(13) 



where T> is the temperature of the uniform matter. 

For finite thermal conductivity, n ^ 0, from Eqs. (fTTj) . (fT2!) we obtain 



Ssn = 6nr 



n[l + fi:p^/(cy7)] 



unpaids /dn) 



7SCy 



The increment, ^{p), is determined by Eq. ( ITOl) . where we put 



(14) 



(15) 



Thus we obtain 



7 



-P 



utp H h cp + 

mn 



ul-ul 



1 + Kp^jicyy) 



(16) 



This equation differs from that has been derived in Ref. [22] by presence of 
extra surface tension term (oc c). 

Eq. ( IT6l) has three solutions. Expanding the solutions for small momenta (long- 
wave limit) we find 



7i,2 = ±iUsP + 



K I Un 



Cy \U\ 



dri + C 



mn 



p" 



73^ 



2 2 
KU^P 

ujcv 



1 - 






C + 



KjUirp 

ulcl 



{drj + Qku^\ 2 



mncyuj 



p 



(17) 
fl8) 



The solutions 71,2 correspond to the sound mode in the long wavelength limit. 
The third solution, 73, is responsible for the thermal transport mode. Below the 
ITS line (and above the AS line) -u^ < 0, m| > 0, solutions 71^2 correspond to 
an oscillation and damping, whereas 73 describes an unstable growing solution. 
Below the AS line, since uf < and m^ < 0, the modes exchange their roles: 
the sound modes become unstable, while the thermal mode is damped. 

We solved the system of non-relativistic hydrodynamical equations. Solutions 
generalizing Eqs. (TT7I) and (ITS]) to the relativistic case are presented in Ap- 
pendix B. 

Further we will consider two limiting cases of a large and small thermal con- 
ductivity, K. 



3.2 Limit of a large thermal conductivity 

For i^p^ /{cv\l\) ^ 1 from Eq. (TT6!) we arrive at 

(19) 



2 2 

1 = -P 



2 , i.dr] + Ql , 2 

Urp H h cp 

mn 



that we have derived in Refs. [T3fT^ . Note that dependence on k disappeared 
in this limiting case. The solutions are 




7 = -p2 / ^ ^^4 _ 1 _ p2y2^. (20) 



Here /3 = c{mnY / {drj+C^Y is the dimensionless parameter. For c\ut\/ 1^ ^ p ^ 



\ut\/\Jc\1 + 1//3| we obtain 7^ = -u^p^, and u^ = -\crV^{T) = 9(5r/o/(4m) 
for the particular case of the Van der Waals fluid in the vicinity of the critical 
point (ricryTcr) (taken as the reference point, cf. Eq. ([6])). 

An instability arises within the ITS region, where u"^ < 0. The most rapidly 
growing mode corresponds to 7 ~ 7^, and p ~ Pm'- 



^"^ (2v^-M)(Jr/ + C)' 



^^ (2v^+l)c- 



2 _ \ "U't 

m 



10 



Using these expressions we rewrite the condition np^ /{cy^) ^ 1 as 

K 3> Cy\fc. 

Then also we have 



nnpl^{ds/dn)T Kn{ds/dn)T 

ImSCv Sy/cCv 

1 + Kpl,/{cvlm) = 1 + ^V' (22) 

Cvy/c 



and following ([T 

(iTo = ario r— — . (23) 

Cyn [1 + K/[y/cCv)\ 



We see that an assumption on the spatial uniformity of the system fails right 
after the ITS region is reached. An aerosol (mist) of bubbles and droplets is 
formed for a typical time t^er ~ l/7m- For the Van der Waals equation of state 
using Eq. (155]) of Appendix A we find that STq/Stiq > 0, i.e. the temperature is 
larger in denser regions. However the amplitude of the temperature modulation 
is rather small for k, ^ Cyy/c. 

In the limiting case, one can neglect at all small deviations of 6Tq and put 
6Tq = in Eq. (I23p . Formally this limit corresponds to setting k — )■ oo in 
Eq. (123 p . In order to avoid misunderstanding we stress that we continue to 
consider only slow motions staying in the framework of the validity of the 
non-relativistic hydrodynamics [28]. 

As before, we assume that at a slow expansion with s{t) ~ const the ini- 
tially uniform system reaches ITS region (see corresponding curves in Fig. 2). 
Solution of Eq. (jlj) for /t — > oo corresponds to T(r) = const, and to the approx- 
imate conservation of the entropy (provided the velocity u or the viscosities 
are small), as it follows from the l.h.s of Eq. (jl]). We may also see it from Eq. 
flTTl) . For that let us divide this equation on n. Then we can put zero the l.h.s. 
of thus obtained equation (since k — )■ oo). Any solution T(r) = const fulfils 
this equation. Using Eq. flT^ we arrive at 6T = (not dependent on r and t). 
The latter approximation was used in Refs. |13fl4j . where we solved Eq. (TTO!) 
for T =const and derived Eq. (TT9|) and Eq. (120|) . On the other hand, from Eq. 
(fT2!) at (5T = we find 6so = {ds/dn)\TSnQ, that follows from Eq. (TT^ as well. 
Using that 6T = in Eq. (IT^ we again arrive at Eq. (IT^ and Eq. (12UI) . 



11 



3.3 Limit of small thermal conductivity 



For the sake of simplicity, here we consider the case, when both shear and bulk 
viscosities are zero, but thermal conductivity is small but non-zero. Dropping 
7^ term in the l.h.s. of Eq. (fT6|) (we justify this approximation below) we 
immediately find the solution 

,, = - '"fp';+<("t (24) 



ujcv 1 + cp'^/u 



Expanding this value in small p we see that it, indeed, corresponds to the 
solution 73 given by Eq. ([H 



The most rapidly growing mode corresponds to 




(25) 
Slightly below the ITS line (for -m^ < 1), 



4 

p^^-4/(2c), 73^ = 73(p^) ^ ^^. (26) 

Now we are able to check the validity of our solution in this case. Comparing 
different terms in Eq. fITBl) we see that it is legitimate to drop the term 7^ in 
the l.h.s. of equation, provided 



K < cvVc—l-. (27) 

Comparing Eq. ( |26|) with Eq. ( 12T|) for /3 — )• 00 (zero viscosity) we see that at 
Cyy/c ^ K <^ Cv\/cuy\u^\, when both solutions 7^ and 73m are valid, the 
following inequality is satisfied 7^ ^ 73m- 

Thus we found that in all considered above cases the instability occurs at the 
crossing of the ITS line. However, for k ^ Cvy/c the most rapidly growing 
mode corresponds to 7™ and in the opposite limit n <^ cyy/c, to 73^- 



12 



3-4 Limit of k, = 

In this specific limit there are significant differences compared to the case 
of any non-zero n. Indeed, for k = there is no solution with 73 7^ 0. Be- 
low the ITS line and above the AS line the thermal mode, that drives the 
system towards equilibrium for non-zero thermal conductivity, does not exist 
for K = 0. Thus the evolution in the spinodal region is entirely governed by 
approximately adiabatic sound excitations. 

From Eq. (ITT]) we find 

6so = 6nQs/n, (28) 

and from Eqs. ([12]) and ([28]) 

Cyn 

Thereby the temperature is modulated, similar to the entropy density. For the 
Van der Waals equation of state using Eq. fl38]) of Appendix A we find that 
STq/Stiq > 0. The amphtude of the temperature modulation is here larger 
than in case of large n, see Eq. ( l23l) . 

The increment 7 is given by Eq. (TT9]) with Ug replaced by ut'- 

(30) 



2 2 

7 = -P 



2 , (c^^ + C)7 , 2 

M; H h CP 

mn 



Thus, contrary to the case k 7^ 0, instability appears, when the system 
trajectory crosses the AS line rather than the ITS line. The value m| = 
—^p,maxv'^{s) = 156 sTp^rnax/ {512m) for the particular case of the Van der 
Waals fiuid in the vicinity of the point {np^max, Tp^max) (taken as the reference 
point, cf. Eq. ([8])). This result holds also in case of the ideal hydrodynamics, 
where in addition to k = 0, the viscosity coefficients {rj, () are zero. 

Note that for any small but finite value k the solution 73 results in instability 
already for -u^ < (i.e. below the ITS line rather than below the AS line). 
Since in reality always k 7^ 0, we conclude that the instability condition m? < 0, 
as one usually derives it within the framework of ideal hydrodynamics, should 
be replaced by the condition m^ < 0. 
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Phase transition from metastable state 



During its evolution a dynamical system crosses a metastable state region 
(binodal) before entering to the region of the instability (spinodal), see Figs. 
1 and 2. If a system evolves rather slowly a phase transition occurs already in 
a metastable region. Therefore let us consider the system in the metastable 
state. Small seeds of the stable phase, being produced in fluctuations, converge, 
while seeds of overcritical sizes start to nucleate, see Ref [31] . In nuclear 
physics these processes have been considered in many works, e.g. see Ref. [2] 
and references therein. 

In Refs. [T^pi] we described the seeds dynamics using equations of non-ideal 
hydrodynamics. Here we continue this description. An analytical description of 
metastable states is possible to perform in the limiting cases of infinitely large 
thermal conductivity (when one can put T =const, see Refs. [T3fT^ ) and of 
zero thermal conductivity. Therefore, below we study the limiting cases at first 
and then we perform estimates for an arbitrary value of thermal conductivity, 

K. 



4.1 Limit o/k — )■ 00 



Let us consider the situation, when at very slow expansion with s{t) ^ const 
spatially uniform spherical system of a large radius Ry achieves either the OL- 
state or the SV- state (see corresponding curves in Fig. 1). During subsequent 
slow expansion, after a while a bubble of an overcritical size R > Rcr of 
stable gas phase, or respectively a droplet of liquid phase may appear in a 
fiuctuation and begin to grow. In case \n — ncr\/ncr ^ 1, i.e. in the vicinity of 
the critical point {ricr^Tcr), corresponding solutions (ITU]1 describing dynamics 
of the density in the fiuctuation were searched in Ref. [13] in the form : 



5n{t,r)-' ""^^^ 



m 



±th -^^^ + 



/ 2A,,t;3(T) 



+ {^n)cor, (31) 



where / = J2c/{Xcrv'^{T)), the upper sign corresponds to the growth of the 
bubble and the lower one, to the growth of the droplet, and the solution 
is valid for \e/{XcrV^{T))\ -C 1, i.e. for 7^ = const -C 1 and for the initial 
bubble/droplet size -R„(0) > RcriT). In Refs. [T^fTl] we neglected a small 
correction {5n)cor in Eq- (I^Tj) . With this small correction one is able to recover 
exact baryon number conservation, see below. One can easily see that Eq. (1311) 
fulfils Eq. ( ITOj) provided the temperature, T, does not depend on time and 
space. As we argued above, this result follows from Eq. ( ITT]) , if one formally 
sets K — )■ 00. 
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In case of expanding system we also assume that typical time for the formation 
and growth of a fluctuation of our interest {tform + tp) is much smaller than 
the typical expansion time texp- 

Substituting Eq. (ETj) to Eq. ( ITU]) and considering r in the vicinity of the 
bubble/droplet boundary we obtain equation describing growth of the bub- 
ble/droplet size [TMi] : 



/3to d Rn _ 3e ^ 2/ ^ to dRn , . 



where to = '^idrjcr + C,cr) /[Xcrv'^{T)'mncr] is the time scale. It follows from the 
solution of this equation that /?„• = 4/Acr'y^(T)/(3e). First the bubble/droplet 
size Rn{t) > Rcr grows with an acceleration and then it reaches a steady grow 
regime with a constant velocity Uas = ^ ^i{T)t ^ l(\Tcr — T|^/^. 

In the bubble/droplet interior 5n ~ ^v(T)/m, and this region grows with 
time. Thus {Sn)cor — v{T)R^{t) /{mR^y) produces compensation of the baryon 
number. Note that this correction is very small for Rn{t) ^ Ry and does not 
affect solution at i? ~ Rn{t)- 

Substituting Eq. f l3T]) to Eq. flT2|) for T = const we obtain 




±ti/ - ^w + 



I 2XcrV^{T) 



+ {6n),or . (33) 



Thus the temperature stays constant in the bubble/droplet interior and exte- 
rior and the entropy and the density differ in the interior and exterior regions. 
In spite of that, approximation of a quasi-adiabatic expansion of the system 
as a whole might be used even, when the system reaches metastable region, 
provided the gas of bubbles/droplets is rare, or the value v{T) is small. 



4-2 Effect of finite thermal conductivity 



Let us determine, when in the case k 7^ we are still able to exploit solutions 
found above for k, — )■ 00. Estimating from Eq. flT2|) 6T ~ 6sTcr/cv,cr (for 
\6s\ ~ \{ds/dn)TSn\) with 6n from Eq. fl3T|) and setting this estimate in Eq. 
(ITTl) we find the typical time for the heat transport t^ ~ R?cv,cr/ i^cr- For 
R{t) < Rfog = K'cr / iuasCv,cr) ouc has t„ > ty and the droplet /bubble dynamics 
is determined by Eq. (!32|) for the change of the density, cf . p^p^ . At k — > 00 
for all times the droplet /bubble dynamics is determined by Eq. (132|) . At k 7^ 
for R{t) > -Rfog, i.e. for t„ < ty, the dynamics is controlled by the heat 
transport. 
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4-3 Limit of k, = 



Limit K = is again specific. In case \n — np^maxl/np^max ^ 1, i-e. in the 
vicinity of the point {np^rnax,sp^max), corresponding solutions (TTUj) describing 
dynamics of the density in the fluctuation can be searched in the form fl3T]) 
with the only difference that 6T should be replaced by 6s, Xcr by Xp^max and 
v{T) by v{s). Thus 



6n{t,r) 



v[sj 
m 



±th -^^^ + 

t 



2A 



P,max 



v^{s) 



+ {Sn) 



(34) 



where / = j2c/[Xp^maxv'^{s)], the upper sign corresponds to the growth of the 
bubble and the lower one, to the growth of the droplet, and the solution is 
valid for |e/[Ap^max'i^^(s)]| -C 1, i.e. for 7^ = const ^ 1 and for -R„(0) > Rcr{s). 
Deriving Eq. (^^ we used that s does not depend on t and r. 



Dynamics of Rn (t) is determined by Eq. f l32|) , where one should replace values 
calculated at {ricT^Tcr) to the corresponding values at {rip^rnax^sp^rnax)- 



From Eq. f ITT]) we obtain 

5s = = {Ssnp^rnax " Sp^jnax 



on)/np^^^ 



(35) 



with 5n given by Eq. flM|) . Knowing 5s and 5n we find 5T from Eq. flT2|) . Thus 
not only density, n, but also the entropy density, s, and the temperature, T, 
change in the surface layer and get different values inside and outside the 
bubble/droplet. Contrary, the value s remains constant. 



5 Dynamical effects of the expansion 



Let \texp{n)]s is the typical time for the expansion of the system (at approxi- 
mately constant entropy) prepared in heavy ion collision and the system tra- 
jectory passes through the ITS region. In case, if taer > [texp{n)]s the aerosol 
has no time to develop in the corresponding density interval. Then the system 
(being covered by fine ripples) continues a quasi-uniform expansion traversing 
the ITS region. Further it may cross the AS region (note that it also can be 
that the expanding system freezes out before the trajectory crossed the AS 
line). As is seen from Fig. 2, the turning line can be further reached, where 
P{nt) = 0. For s > Sa the turning points are lying inside the AS region. For 
St2 < S < Sti the turning points are outside the AS region, but within the 
ITS region. In the vicinity of the turning point, the expansion time t(,xp{nt) 
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can be significantly increased yielding taer < texp{nt)- Note that, as follows 
from the analysis with the realistic equations of state the turning points can 
be reached only for low values of the entropy s. Such low values are reachable 
in case of low energy heavy ion collisions, where there arises nGL transition. 
In case of violent heavy-ion collisions at conditions of the first-order hQGP 
phase transition values of s are much larger and P is always positive. 



6 Concluding remarks 



In this work we found various solutions describing dynamics of the density, 
temperature, and entropy density in the spinodal and metastable regions at 
first-order phase transitions. We tried to answer the question how one should 
treat the fact that critical temperatures at the isothermal and adiabatic lines 
are different. In particular, we demonstrated that correct region of instability 
is not appropriately found within standard ideal hydrodynamical calculations 
(where instability manifests itself at the crossing of the adiabatic line). This 
might be an important finding since ideal hydrodynamics is often used in 
actual simulations of heavy ion collisions. Since in reality k 7^ 0, spinodal in- 
stabilities take place, when the system phase trajectory crosses the isothermal 
spinodal line in the course of the heavy ion collision. The value of the critical 
temperature, Tcr-, at the pressure isotherm [P{n) for fixed T] is significantly 
larger than the value of the maximum Tp^^ax at the pressure adiabat [P{n) 
for fixed s]. This means that according to our findings the prospects of observ- 
ing signatures of spinodal instabilities in experiments with heavy ions become 
more promising. 
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Appendix A. Van der Waals equation of state 

All thermodynamical quantities can be found from expression for the free 
energy 

F[n, T] = Nfo [t [in(A^n) - l] - T ln(l - bn) - an} , (36) 
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where A = ^271 /{mT) is the thermal wave-length, N is the baryon number. 
Parameter a governs the strength of a mean field attraction and h controls a 
short-range repulsion, /q is an extra constant, which can be fitted to obtain 
more realistic values of thermodynamic quantities for the specific phase tran- 
sition (e.g., for nGL or hQGP). From Eq. (I36ll we find the pressure in {n,T) 
variables: 



P[n, T] = /o nT/{l - hn) - n^a 



(37) 



and the entropy per baryon 

s[n, T] = /o [- ln(AV6) + ln(l/(6n) - 1) + 5/21 + C^. 



(3J 



Since entropy in classical systems is defined up to a constant, an arbitrary 
constant Ci is added. 

From Eqs. (1371) and (l38|l we also find the energy 

E[n, s] = Nfo{3Bx'^/'^ - an), (39) 



the pressure 

P[n,§] = /o (25x^/3 -an^). 



(40) 



and the temperature 
T[n,5] =25x2/^ 



(41) 



expressed in terms of {n, s) variables. Here B = (7r/m)exp[2s/(3/o) — 5/3 + Ci], 
X = n/{l~bn). In numerical calculations, which results are illustrated by Figs. 
1,2, being performed just with a demonstration purpose we put /o = 1 and 
Ci = 0, as for the case of the standard Van der Waals equation of state. 

Appendix B. Spinodal instability in relativistic hydrodynamics 

In this Appendix we generalize results obtained in Sect. [3] to the relativistic 
case. We use the Landau approach [32] • It is well-known that equations of 
relativistic first-order hydrodynamics are parabolic and formally suffer of the 
causality problem. However the causality problem may arise only, when one 
tries to describe the fiow with a fast velocity and for short wavelength phe- 
nomena. To avoid this problem following Ref. [33] we assume that the velocity 
of the collective motion is much slower than the mean particle velocity in the 



system (the latter can be relativistic) and we are interested in description of 
a long wavelength phenomena. 

Let the equilibrium state be the fluid rest frame. Then 5u^ = {0,6u). As in 
non-relativistic case we linearize the hydrodynamical equations with respect to 
the variations Sn, Ss, 6T and Su. With the help of the standard thermodynamic 
relations the Landau equations are reduced to [33 j : 



TlT f 11) 

dtSn + niV5vi) - k^-V^ 6P - cnV^Sn 6T 

w^ \ T 

wdt{V6u) - V^ [(C + rfr7)(V5u) - 6P + cnV'^dn 

dt5s - -dt6n + -V^ (sP - cnV^6n - -6T 
n w \ T 

where w = e + P and s are the enthalpy density and the entropy density, 
respectively. Here we additionally included surface tension effects. 

The pressure and the entropy density variations are expressed through the 
density and the temperature variations (see e.g. Eq. flTSj) for the pressure 
variation). Substituting Eqs. flT^ in Eqs. (I42ti44p we flnd the dispersion law. 
In the long wavelength limit (up to p^-terms) we obtain 



j^o. 


(42) 


= 0, 


(43) 


0, 


(44) 



Yii = ±iCsP + 




(C + dv) 



^ '''^■{2a - !)«?-, (45) 



w ^ '2 



iT = -^^P'- (46) 

upv 

Here c| = {dP/de)s = {n/w){dP/dn)s = {mn/w)uj, a = cy/cp = u^uj'^. 
Similar to the non-relativistic case, the solutions jf'l describe sound excita- 
tions in the long wavelength limit, while 73*^' corresponds to the thermal mode. 
The last term in expression for 'yl'^l is a small correction oc O (T/m) for T <^m 
of our interest. Dropping it we arrive at Eqs. ( TT71) . ( IT8|) . which we have derived 
in the non-relativistic limit in the paper body. Thus we come to precisely the 
same conclusions about instability region, as in the non-relativistic case. 
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